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We examine the single-particle excitation spectrum in the one-dimensional Hubbard-Holstein 
model at half-filling by performing the dynamical density matrix renormalization group (DDMRG) 
calculation. The DDMRG results are interpreted as superposition of spectra for a spinless carrier 
dressed with phonons. The superposition is a consequence of robustness of the spin-charge separation 
against electron-phonon coupling. The separation is in contrast to the coupling between phonon 
and spin degrees of freedom in two-dimensional systems. We discuss implication of the results of 
the recent angle-resolved photoemission spectroscopy measurements on SrCu02. 
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The interplay between electron correlation and 
electron-phonon coupling is one of the hot topics in 
the field of strongly correlated electron systems such as 
high- Tc cuprates. Particularly in an electron-removal 
process from the Mott insulators, the electron-phonon 
coupling occurs due to charge imbalance around a cre- 
ated hole. Thus, the angle-resolved photoemission spec- 
troscopy (ARPES) is a direct tool, and observes a 
quasiparticle dressed with phonons 1]. For the two- 
dimensional (2D) insulating cuprates, the ARPES ex- 
periments have revealed a broad low-energy peak that 
is interpleted as a result of disappearance of the quasi- 
particle weight due to the coupling H, Q . 

The quasiparticle in the 2D systems is not only dressed 
with phonon cloud, but also dressed with antiferromag- 
netic (AF) spin fluctuation. Thus, the effect of phonon 
on the spectrum is affected by the spin configuration of 
the background. However, in one- dimensional (ID) cor- 
related electron systems, a photohole created by ARPES 
decays into spinon and holon due to the spin-charge sep- 
aration 0. Therefore, the effect of phonon on the spec- 
trum strongly depends on whether the spin-charge sepa- 
ration is robust against the electron-phonon coupling. In 
this Letter, we examine the effect of phonon on the single- 
particle excitation spectrum in ID Mott insulators. 

Theoretically, the Hubbard-Holstein model is a basic 
model to study the interplay between electron correlation 
and electron-phonon coupling in cuprates 0, Q. How- 
ever, we have only limited information on the sin gle- 
particle excitation spectrum in this model 0, 0, IS IS 
im 01 . This is because it is hard to treat the infi- 
nite phononic degrees of freedom and electron correlation 
on an equal footing. In order to overcome the difficulty, 
we apply the dynamical density matrix renormalization 
group (DDMRG) method to the calculation of the spec- 
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tra in the half-filled systems. 

We find the following four characteristic features in 
the single-particle excitation spectrum: (i) a dip at high- 
binding energy side of the spinon branch, (ii) broad holon 
branch, (iii) decrease of the spectral weight of the spinon 
branch, and (iv) slight enhancement of the weight at low- 
binding energy side of the spinon branch. In the light of 
the spin-charge separation, these results are reproduced 
well by performing superposition of the spectra for a spin- 
less carrier dressed with phonons. We examine the differ- 
ence between the phonon effects on ID and 2D systems. 
Implication of the results for the recent ARPES measure- 
ments on a ID Mott insulator SrCu02 is discussed. 

The Hubbard-Holstein Hamiltonian in ID is defined by 

H = -iI](cUc.:+i,^+H.c.) 

i.(7 




+uoJ2b+h~gY,ibt + b^in^-l), (1) 
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where {ci,a) is the creation (annihilation) operator 
for an electron with spin a at site i, bf {bi) is the cre- 
ation (annihilation) operator for an Einstein phonon at 
site i. Hi = Ui^f -I- rii^i, rii^a = cJ^Ci^o-, t is the hopping 
integral, U is the on-site Coulomb repulsion, is the 
phonon frequency, and g is the electron-phonon coupling 
constant. We take U to he U = 8t, which is an appro- 
priate value for the cuprates. Here, we briefly mention 
how to construct the model Hamiltonian in the case of 
the cuprates. The lower Hubbard band of Eq. ^ corre- 
sponds to the Zhang-Rice singlet which is derived from 
a three-band model for Cu 3(1^2 _y2 and O 2p orbitals 
with lattice distortion. The distortion leads to change of 
the hopping integral of an electron between neighboring 
Cu and O orbitals. Due to the modulation of the hop- 
ping integral, the diagonal (Holstcin) term dominates in 
the effective Hamiltonian for the Zhang-Rice singlet in 
comparison with off-diagonal terms |^ Q ■ 
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We examine the single-particle excitation spectrum at 
zero temperature defined by 



A{k,uj) = lm(0 



^''''^ En - UJ - H 



Z7 



-Cfc,t 



0), (2) 



where Ck^-\ is the momentum representation of the elec- 
tron operator Cij, |0) denotes the ground state with en- 
ergy Eq, and 7 is a small positive number, which is taken 
to be 7 = O.lt in the present calculation. 

Here, A{k,uj) is calculated by the finite-system 
DDMRG algorithm jM, El 111 El El . The sys- 

tem size L is taken to be 20 lattice sites. The DMRG 
bases are truncated up to m = 400 states from the den- 
sity matrix for a mixed state of the ground state |0), the 
final state after the one electron-removal process Ck,^ |0), 
and two correction vectors {Eq — lu — H + i^)~^Ck,i |0) 
for w = ijji,uj2, and UJ2 — uji — 27. It is technically use- 
ful to note that the convergent results are obtained when 
A{k,uj) and A{k,uj + 27) are smoothly connected for a 
given m. 

keep numerical precision and reduce 
, an open boundary condition is taken 
—trii at the edges. In such a system, 
k is defined by k = mr/{L + 1) with 
The momentum representation of is 

lilEi 



In order to 
boudary effects 
with potentials 
the momentum 
n = 1,2, ...,L. 
given by Cfe,| = ^2/(i + 1) Y,i sm{kl)ci^i 



It is 



noted that the spectrum for L = 20 presented here is sim- 
ilar to that for L = 120 without the edge potential Ell- 
in the DMRG calculation, electronic and phononic de- 
grees of freedom at each site are devided into two differ- 
ent pseudo-sites. Furthermore, the phonon site is decom- 
posed into a set of N hard-core bosons, where the phonon 
states at each site are truncated up to M = 2^. The ex- 
act transformation between the ordinary boson operator 
and the set of the hard-core boson operators has been 
discussed in ref [2H . After the transformation, the high- 
est order term of the bosons is renormalized at first in 
some processes. The processes worsen numerical preci- 
sion. Thus, we set up superblocks so that the final sweep 
process renormalizes the bosons in turn. Here, the max- 
imum number of A'^ which we take is 4. The superblock 
is finally composed of (1 -I- N)L = 100 pseudo-sites. 

Figure 1 shows A(k, u) with and without the electron- 
phonon couping g. The phonon frequency is taken to 
be loq = Q.bt > 7 = O.lt in order to see the effect of 
phonon clearly. However, the frequency for the cuprates 
is Wo ~ 7. Thus, we have confirmed that our conclusion 
does not change for ojq = Q.2t. The origin of energy is 
located at the center of the Mott gap. In Fig. 1(a), two 
branches disperse, merging toward k — > 7r/2. These en- 
ergy positions are us/t = —3.13 and —4.90 at fc = 7r/21, 
which is the minimum momentum in the calculation. The 
branch located in low (high) binding energy side is de- 
duced to be the spinon (holon) branch. Here, fine struc- 
tures inside the branches come from the finite size effect. 
It is noted that the band width of the holon branch from 
k = 7r/21 to k = 107r/21 depends on the open bound- 
ary condition and U. In Fig. 1(b), we find that the 




FIG. 1: A{k,iu) for the ID Holstein-Hubbard model at half- 
filling. The momentum is taken from n/21 to ll7r/21. 
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FIG. 2: (a) A{k 0, ui) for various uiq values. For luq ^ 0, g 
is taken to be 0.5t. (b) A{k — > 0, u) with loo = 0.5t for various 
g values. For (a) and (b), the spectra with g are shifted 
so that the energy of the spinon branch is taken to be equal. 



holon branch becomes broad due to the electron-phonon 
couping, while the spinon branch is sharp. We also find 
a 'peak-dip-hump' structure at the high-binding energy 
side of the spinon branch. The dip disperses like the 
spinon branch. As shown in Fig. 2(a), the energy of the 
dip (hump) position decreases with increasing luq. The 
energy difference between the dip (hump) and the spinon 
branch is estimated to be luq, which means that the peak- 
dip- hump structure is due to the phonon effect. There- 
fore, the orign of the peak-dip-hump structure provides 
information on the effect of phonon. 

As g goes from 0.5t to larger values with keeping 
luq = 0.5t, the spinon branch starts to broaden [oj. In 
Fig. 2(b), we show A{k 0,w) for various g values. 
The weight of the spinon branch decreases linearly as a 
function of g for g > luq. For g = t, the holon branch is 
completely smeared out. The peak-dip- hump structure 
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FIG. 3: (a) Shadded area: schematic view of the dispersion 
in a ID Mott insulator; dashed (solid) line: one holon dis- 
persion for g — {g ^ 0). A dotted line guides the eye. (b) 
A-osik = 0,w); dashed (solid) line: g — {g = ujo ~ 0.5t). 
Here, the broadening 7 is assumed, (c) Ah{p — tt, e) for the 
ID Holstein model, i.e. a model for a spinless carrier with 
Einstein phonons {g = loq — 0.5t). A dashed line is a cosine 
band with the width At. 



observed in Fig. 1(b) develops into multipeaks whose 
positions are uj/t — —3.13, —3.61 ^ —3.13 — luq, and 
—4.09 ^ —3.13 — 2luo, respectively. Furthermore, a shght 
enhancement of the low-binding energy side of the spinon 
branch is seen. The energy difference between the en- 
hancement and the band top {k = 7r/2) for g = is 
estimated to be ojq. Therefore, the three main peaks and 
the enhancement are also due to the phonon effect. 

As mentioned in the previous two paragraphs, the fol- 
lowing phonon effects appear on A{k,uj): (i) the peak- 
dip-hump structure, (ii) the broad holon branch, (iii) the 
decrease of the spectral weight of the spinon branch, and 
(iv) the enhancement of the low-binding energy side of 
the spinon branch. In particular, the characteristic en- 
ergy scale of (i) and (iv) is ajQ. In the light of the spin- 
charge separation, they can be interpreted as follows 

First, we consider the origins of (i) and (ii) shown 
in Fig. 1(b) for g = 0.5t = luq. Let us start with 
the dispersion in the Hubbard model illustrated in Fig. 
3(a). According to the Bethe anzats solution, the disper- 
sion is constructed by supersposition of a set of holon 
dis per sions forming a cosine band with the width of 
At [221 . The superposition is a consequence of the spin- 
charge separation, because each of the holon dispersions 



is characterized by one spinon momentum. Therefore, 
an effective model of A{k,uj), Aeg{k,(jj), is constructed 
by putting the spectral weight for a spinless fermion, 
Aii{p,e) — 5{e — 2tcosp), on each of the holon disper- 
sions. Since a top of the consinc band is running along 
the spinon dispersion es('7+7r/2) = — (7rJ/2)| sin((7+7r/2)| 
for — 7r/2 <q< tt/2, AcB{k,uj) is defined by 

7r/2 

A^s{k,uj)= ^/i(fc-9,'^ + 2i + e,(g+|)),(3) 

except for constant energy shift. A dashed line in Fig. 
3(b) shows Acff{k = 0,lu). The singularity of the spinon 
branch appears aX lu = —ttJ/2 ~ —2Tit^/U = — 0.785i. 
The singularity comes from the flatness of the spinon dis- 
persion near k — Q. The lineshape is consistent with the 
DDMRG data in Fig. 1(a) except for the singularity of 
the holon branch. The singularity of the holon branch is 
recovered from the phase string effect [2^. The consis- 
tency indicates that Eq. Q is appropriate for the spectral 
weight in spin-charge separated systems. 

Let us introduce the electron-phonon coupling, and 
take g to be 0.5t. Due to the spin-charge separation, 
each of the holons couples with phonons independently. 
Namely, Ah{p,e) is given by the spectra for a spinless 
carrier dressed with Einstein phonons. Figure 3(c) shows 
Ah {p — TT, e) that splits into low- lying peaks and an in- 
coherent part 25^ 2^ 2J} ■ The split occurs at the 
anticrossing point e ~ 1.5t that is away from the top of 
the band (p = tt) by wo. At e ~ 1.5t, we find a tiny spec- 
tral weight with a flat dispersion coming from the phonon 
branch. In Fig. 3(a), the split of one holon dispersion 
is illustrated. Aofr(A: = 0,0;) is then given by a solid line 
in Fig. 3(b). A peak-dip- hump structure appears. The 
spectral weight lost by the dip is transfered to high en- 
ergy region. The spinon branch, the dip and the broad 
holon branch originate from the low-lying peak, the anti- 
crossing and the incoherent part of Ah {p, e) , respectively. 
These features obtained by Eq. Q are consistent with the 
DDMRG data, and thus the spin-charge separation is ro- 
bust. 

Next, we consider the origin of (iii) shown in Fig. 2(b) 
for g > loq. As mentioned in the previous paragraph, the 
spinon branch can be expressed by the superposition of 
the low-lying peak in Fig. 3(c). It has been shown that 
the spectral wei ght of the low-lying peak decreases with 
increasing g ^HHM 113 ■ Therefore, the weight of the 
spinon branch decreases with increasing g. 

Finally, the origin of (iv) is considered. In Fig. 4, 
we see that the origin comes from the phonon branch. 
Even for the coupling g — t > luq, the spectral weight of 
the phonon branch is very weak at fc [23. 125I l2a. l27l | . 
This is the reason for the tiny weight seen in the DDMRG 
data. It is noted that the tiny weight is observed when a 
condition ujq < tt J/2 is satisfied. 

All of our interpretations for (i)-(iv) show robustness 
of the spin-charge separation. The robustness leads to 
the difference of the effect of phonon between ID and 
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FIG. 4: (a) Shadded area; schematic view of the dispersion in 
a ID Mott insulator; dashed (sohd) hne: one holon dispersion 
for g = {g ^ 0). A dotted line guides the eye. (b) A{k 
0, ui) for g — t and ujo ~ 0.5t. 



2D systems in the presence of electron correlation. In 
order to see the difference, we introduce a dimensionless 
parameter A — g^juioW with the noninteracting band 
width W . There is a characteristic A value, A* |23|. For 
A < A*, the lowest energy excitation is weakly dressed 
with phonons. As A approaches to A*, the excitation 
loses its weight rapidly. Then, the dominant low-energy 
excitation moves to a heavily dressed polaron. In the ID 
Holstein model. A* is estimated to be A* 1 [29l Isoj. 
In the ID Hubbard-Holstein model, the spinon branch 
loses its weight for A* ~ 1 (g* ~ \At) as estimated 
from Fig. 2(b) 0. This is because the spinon branch 
can be expressed by the superposition of the spectra for 
the low-energy excitation of the Holstein model. In 2D, 
on the other hand. A* is close to 0.2 in the i- J-Holstein 
model with J = 0.3t, while A* > 0.6 in the Holstein 
model 0, 0, 1^ . This means that the AF correlation 
helps the formation of the lattice polaron [s^ls^- ls^ . As 
J increases, the energy gain is caused by the AF correla- 
tion. Then, the coherent motion of a hole is suppressed. 



The suppression helps the formation of the lattice po- 
laron. 

Finally, let us discuss the ARPES data for a ID Mott 
insulator SrCu02 in the hght of the DDMRG data. In 
this compound, high-energy ARPES experiments have 
been done, where the spinon and holon branches were 
observed j5]. Near the P point, the intensity of the 
holon branch is smaller than that of the spinon branch. 
In addition, these branches do not exhibit singularities 
predicted by the Hubbard model [S^. In Fig. 2(b), we 
show the DDMRG data for A = 0.25 (g = 0.707t and 
cjQ = 0.5t) that is appropriate for ID cuprates 7j. The 
DDMRG data is consistent with the ARPES data except 
for the dip. The dip may be covered by the phonon dis- 
persion. Finite temperature effects in the presence of g 
might be the reason for the covering. 

In summary, we have examined the single-particle ex- 
citation spectrum in the ID Hubbard-Holstein model at 
half-filling by performing the DDMRG calculations. We 
found the following four characteristic features: (i) peak- 
dip-hump structure, (ii) broad holon branch, (iii) the de- 
crease of the spectral weight of the spinon branch, and 
(iv) the enhancement of the weight at low-binding en- 
ergy side of the spinon branch. In the light of the spin- 
charge separation, these results are reproduced well by 
performing superposition of the spectra for a spinless 
carrier dressed with phonons. In ID systems, the spin- 
charge separation is robust against the electron-phonon 
coupling, while in 2D systems, the AF spin correlation 
favors the formation of the lattice polaron. The ARPES 
data for SrCu02 was consistent with our DDMRG data. 
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